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Abstract 



I develop n-site cluster approximations for a stochastic sandpile in one di- 
mension. A height restriction is imposed to limit the number of states: each 
site can harbor at most two particles (height zi <2). (This yields a consider- 
able simplification over the unrestricted case, in which the number of states 
per site is unbounded.) On the basis of results for n < 11 sites, I estimate 
the critical particle density as C,c = 0.930(1), in good agreement with simula- 
tions. A coherent anomaly analysis yields estimates for the order parameter 
exponent [/? = 0.41(1)] and the relaxation time exponent {u\\ ~ 2.5). 
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I. INTRODUCTION 



Sandpile models are the prime example of self-organized criticality (SOC) in which 
a system with an absorbing-state phase transition is forced to its critical point , leading 
to scale-invariance in the apparent absence of parameters [Q. The absorbing-state phase 
transition, which depends, as is usual, on the fine-tuning of one or more control parameters, 
is evident in sandpiles with a fixed number of particles [^fzHTOl] , models that have come to be 
called fixed-energy sandpiles (FES). While most studies of sandpiles have probed the driven 
case P,p!^, there is great interest in understanding the scaling properties of FES models as 
well g^THH- 

Previous studies of FES reveal that they exhibit a phase transition between an absorbing 
and an active state as the particle density ( (which is the temperaturelike control parameter) 
is increased beyond a critical value . Until now, all quantitative results on FES have 

been obtained from simulations. It is therefore of interest to apply theoretical methods to 
such models. One such approach is Suzuki's coherent anomaly method (CAM) for analyzing 
a series of cluster approximations. It has been shown to yield good estimates for critical 
properties both in [0,0 and out of equilibrium |T9|J20[] . In this work I develop n-site 
approximations for a one-dimensional sandpile model, and analyze the results using the 
CAM. This represents the first application of the CAM to a model representative of the 
class of absorbing-state phase transitions in systems with a conserved density pT|-p^. 

In this paper I study a FES with a height restriction. From the theoretical viewpoint, an 
inconvenient feature of sandpile models is the unbounded number of particles that may oc- 
cupy the same site; this complicates attempts to derive cluster approximations. In Manna's 
stochastic sandpile |^,^|, sites with height z > 2 are active. If we restrict the height (or 
number of sand grains per site), to be < 2, the effect on critical properties should be minimal, 
aside from a possible shift in the critical density (c- This expectation was recently verified 
numerically: the restricted-height stochastic sandpile belongs to the same universality class 
as its unrestricted counterpart |^6[. I study the restricted-height model for calculational con- 



venience, as a representative of a broader universality class that includes Manna's stochastic 
sandpile, the conserved lattice gas, and the conserved threshold transfer process ^1 . 



The balance of the paper is organized as follows. The model is defined in Sec. II, followed 
by a discussion of n-site approximations in Sec. III. Numerical results are presented in Sec. 
IV. The CAM analysis is discussed in Sec. V, and in Sec. VI I present a brief summary. 



II. MODEL 



The model is defined on a ring of L sites with periodic boundaries. (The cluster approx- 
imations effectively study the L ^ oo limit.) The configuration is specified by the number 
of particles 2;j = 0, 1, or 2 at each site; sites with Zi = 2 are said to be active, and have a 
toppling rate of unity. The continuous-time (sequential), Markovian dynamics consists of 
a series of toppling events at individual sites. When site i topples, two particles attempt 
move to randomly chosen nearest neighbors j and j' oi i. The new position of each particle 
is accepted if and only if the target site has fewer than two particles. I consider a stochastic 
toppling rule in which the two particles move independently. Any particle attempting to 
move to a site already harboring two particles is sent back to the toppling site. (Thus an 
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attempt to send two particles from site j to site k, with Zk = 1, results m. Zk = 2 and Zj = 1.) 
Transition probabilities are listed in Table I. 

This model, and a closely related one (with a cooperative toppling rule), were studied 
via simulation in Ref. p6[, which showed that the critical exponents (3 and z/_|_ are the same 
as for the unrestricted Manna sandpile. A similar conclusion was reached in Ref. |21[] for a 
two-dimensional restricted-height sandpile (called the conserved threshold transfer process 
in that work). 



III. CLUSTER APPROXIMATIONS 

I have derived dynamic n-site cluster approximations for the one-dimensional restricted 



sandpile model. Such approximations often yield qualitatively correct phase diagrams p7 
The procedure parallels that used by Ferreira and Mendirata to study the one-dimensional 
contact process [^]. The n-site approximation consists of a set of coupled differential equa- 
tions for the probabilities P^'^^ of each ?7,-site configuration, C. (There are 3" such configu- 
rations, but the number of independent probabilities is ~ 3"/2, due to symmetries.) The 
system is assumed homogeneous, so that the P^"^ ^^^^ independent of position. 

Since transitions in a set of n contiguous sites generally depend on sites outside the 
cluster, the n-site probabilities are coupled to those for n+1 and so on, generating an infinite 
hierarchy of equations. The n-site approximation truncates this hierarchy by approximating 
m-site probabilities (for m > n) in terms of n-site conditional probabilitites. In the n-site 



approximation, the joint probabihty for a sequence of n+1 sites is approximated so |£0 



P("+l)(Zi,...,^„+l) ~P(")(z„+l|^„,...,^2)P^"n-^n,...,^l) 
_ P^^\zn+l,...,Z2)P^''\zn,...,Z,) 
P(-1)(Z„,...,Z2) 



(1) 



The equations for one- and two-site approximations are relatively simple to derive, and 
are described in Ref. p6|. I have developed a computational algorithm capable of generating 
the approximation for arbitrary n. Each configuration C = {zn, ■■■,zi) is represented by an 
integer 

n 

m = Y.z,. t-' . (2) 

k=l 

The calculation begins with the generation of all configurations, corresponding to each 
integer from zero (all sites empty) up to the maximum, 3" — 1 (all sites doubly occupied); 
the symmetry (under inversion) of each configuration is determined. If C is not symmetric, 
then it and its mirror image Cr must have the same probability, and only the smaller of C 
and Cr is treated explicitly, reducing the number of variables by roughly half. 

Next, a list of all possible transitions is constructed. Here it is useful to distinguish 
between centra/ transitions (involving a toppling at one of the sites 2,...,n— 1) and boundary 
transitions, in which either site 1 or site n, or one of the peripheral sites (0 or n+1) topples. 
The rate of a central transition C ^ C is the product of a branching probability (for the 
particles to be redistributed in a particular manner, as in Table I), and the intrinsic toppling 
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rate, which is unity. Consider, for example, configuration C = (21120). The transition rates 
associated with a topphng at the second site (counting from the right), are 

w[(21120) ^ (21102)] = ^, 
w[(21120) ^ (21201)] = ^, 

w[(21120) (21210)] = i. 

For each configuration C, the set of allowed transitions to other states, and the associated 
rates, are stored. In the case of a central transition from C to C, the contribution to the 
time-derivative of the probability has the form: 

"^^^^'^ = w\C ^ C']P(C). 



dt 

There is of course a corresponding loss term for P{C): 

dP{C) 



dt 



-w{C C']P{C). 



Thus central transitions contribute to the evolution of the probability distribution precisely 
as in the master equation. 

For boundary transitions, one does not have access to the n + 1-site probabilitites re- 
quired to mount a complete description, and so must resort to the truncation scheme em- 
bodied in Eq. (|lD. For example, the contribution to dP{2, Z2, Zn)/dt due to the transition 

(2, Z2, Zn) —>■ (1, Z2, Zn) is 

ip("+^)(0, 2, Z2, Z^) + ^P("+l)(l, 2, Z2, z^) . 

The p("'+^) are estimated using Eq. (|l]). For boundary transitions one stores not only 
the rate, but the two configurations (aside from the original one, C) whose probabilities 
are needed to evaluate dP{C)/dt. With this information available, one can evaluate the 
derivatives dP{C)/dt for all possible configurations, given the probability distribution. 

The evolution of the probability distribution is found via numerical integration, using a 
fourth-order Runge-Kutta scheme |2^. The integration is halted when a stationary distri- 
bution is attained, that is, when the time-derivatives dP{C)/dt all have an absolute value 
smaller than S (typically, 6 = 10~^^). An interesting technical point concerns the evaluation 
of the n — 1-site marginal distribution. There are evidently two equivalent expressions that 
may be used: 

P^^'-'Hzr,., Z,) = Y P(") 



n—li ■■■■! Zi) / ^ J \Zni Zn—li Zi J 
:n=0 



and 
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20=0 

Numerical stability is greatly enhanced using the mean of the two expressions given above. 

For sizes n > 7, very near the critical point, relaxation to the stationary distribution is 
very slow, and the following procedure proves advantageous. Let D = niax^ \dP{C)/dt\ be 
the largest derivative (in absolute value). The properties of interest (principally, the active- 
site density) are recorded as a function of D, and the integration halted when D < 10~^°. 
Fig. 1 shows the result for the active site density, for ( slightly above The stationary 
value is obtained via extrapolation to = 0, usually via a quadratic fit to the four data 
points for smallest D. (The resulting correction is typically less than 1% of the value at 
D = 10~^°.) I also studied the order parameter relaxation rate 7 = \p/p\ for each n at a 
series of ( values near, but below, Cc- These data are used to estimate the critical exponent 

in Sec. V. 



IV. NUMERICAL RESULTS 

I derived cluster approximations for n < 11 sites, yielding the stationary active-site 
density function of (. The n-site approximation predictions for Pa(C) (^ = 3 to 11) 

are compared against simulation in Fig. 2; the theoretical curves appear to approach the 
simulation result systematically. 

For each n, the active-site density is zero below a certain critical value, (c,n- Since the 
phase transition in the stochastic sandpile is continuous |T§J2^ , one expects the same to be 



true of the cluster approximations. This is indeed the case for n < 4, but for n = 5, 6 and 
7 there is a very small discontinuity in pa (~ 10~^, invisible on the scale of Fig. 2), as we 
decrease (. Since the same procedure is used for all n, the discontinuity is unlikely to be 
artefact of the numerical method. On the other hand, I do not regard the discontinuity as 
physically significant; it appears to represent an unphysical feature of the cluster approxi- 
mations for certain n values, very near the critical density. In the CAM analysis I disregard 
the behavior of pa in the immediate vicinity of (c,n, and instead analyze its properties at 
points somewhat removed from the transition. 

There remains, naturally, the problem of estimating the critical density, (c,n- For each n 
I determine the critical density by fitting the four or five data points nearest the transition, 
where p ~ 10~^ or less. In each case, I plot In versus ln{( — Cc,n); varying (c,n to obtain 
the best power-law fit. The associated slopes vary between 1 (for n = 3) and about 0.25 
(for n = 6), but these, again, are regarded as unimportant details of the approximation 
in question. It is important to stress that, in the cases where the transition is apparently 
discontinuous, the difference between the location of the discontinuity and the extrapolated 
value of (c,n is less than one part in 10^, and that the estimates for and critical exponents 
are insensitive to these tiny differences. The values of (c,n obtained in this manner are listed 
in Table II. 

Using the results for (c,n, I estimate Cc = lini„^oo Cc,n by plotting A„, = Q ~ Cc,n versus 
n in a double- logarithmic plot, varying (c to obtain the best power-law fit. The latter is 
obtained using (c in the range 0.929 - 0.931, yielding (c = 0.930(1), in good agreement with 



the simulation result of 0.92965 I^Gf. The finite-size scaling prediction for the critical point 
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shift is f 251: A 



found in simulations 



(xl/n. I obtain a good fit to the data (see Fig. 3) using z/j_ 
1^.2^), including a correction to scaling term: 



1.66 fas 



A 
oc — 

n 



B 



n 



The numerical data are consistent with the simulation estimate for i^^, but not sufficient to 
furnish an independent estimate of the exponent. 



V. COHERENT ANOMALY ANALYSIS 

A detailed explanation of the CAM procedure is given in Ref. [|17| , |18| ; it may be un- 
derstood on the basis of finite-size scaling The approach here parallels that used by 
Tome and de Oliveira in their study of the Domany-Kinzel model |^0|. To begin, one argues 
that the cluster size n plays the role of an effective system size L as regards scaling proper- 
ties. This is because the n-site approximation effectively cuts off correlations of range > n 
(notwithstanding the fact that cluster approximations nominally treat an infinite system). 
Thus, as noted above, one expects a critical point shift A„ oc n'^^'^^. Finite-size scaling 
theory also yields the relation pa,niCc) oc n'^^"^ oc for the order parameter in a finite 
system, at the (true) critical point. For Q > Cc,n, Pa,n(C) is a smooth function. Thus we are 
led to a scaling hypothesis for the order parameter pO| : 

PaAO = Kf (^) , (3) 

where f{x) is a scaling function with /(O) =0. If we suppose that f{x) oc x^'^'^ for < x < 1, 
then Pa,n{0 = An{C — Cc,nY'^"' , whcrc the amplitude An diverges as n oo: 

An OC A;^^'^^^-^). (4) 

This is the usual CAM relation. On the other hand, the hypothesis that n is equivalent to 
a finite system size leads directly to: 

p,,„(C) ~ A^. (5) 

This expression involves the behavior of the n-site approximation at the critical point C,c 
not Cc,n- It is interesting to note that the hypothesis of an effective system size directly 
proportional to n is not strictly necessary. The scaling relations involving A„ follow from 
the more general hypothesis of an effective system size Le// = Le//(n), for example Le// oc n"^ 
with > 0. 



A. CAM analysis for [5 

As noted above, the n-site approximations for the order parameter pa^n are not all well 
behaved in the vicinity of Cc,n- For this reason, analysis of pa^n at Cc,n will not yield a 
consistent set of well defined amplitudes But since Pa,n{C) is well behaved for ( > (c,n, 
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we can study its scaling at some point intermediate between (c,n and Cc- In particular, the 
scaling hypothesis Eq. (^) implies that if we fix x = {( — (c,n)/'^n, then 

^ = /(x)A^/^-. (6) 

Our strategy is to analyze the order parameter data reasonably near the n-site critical value, 
but away for (c,n itself, where pa,n is singular. A crucial point in this analysis is the postulate 
that the mean-field exponent Pmf = 1, regardless of the behavior of pa,n "in the immediate 
vicinity of (c,n- The motivation for this assumption is, firstly, that Pmf is clearly unity for 
n = 1, 2, or 3; secondly, that a critical exponent such as /? is determined, in mean-field 
theory, by symmetry properties of the order parameter, and hence should not vary with n; 
and thirdly, that Pmf = ^ generically for phase transitions to an absorbing state [^. (The 
basis for this last assertion is that the mean-field equation for the order parameter will have 
the form dp/dt = Ap — Bp^, barring some coincidence or a symmetry that renders A and/or 
B zero |3|.) 

I evaluate dpa,n/dC (numerically, using an interval A = 0.0005), for fixed x = 1/4; the 
results are shown in Fig. 4. Least-squares linear fits to the data for = 8 - 11 yield, via 
Eq. (^, the value (3 = 0.408(6), where the figure in parenthesis denotes the uncertainty. 
Using Cc = 0.929 instead of the best estimate, 0.930, I find (3 = 0.421(5). Thus a reasonable 
estimate for (3 is 0.41(1). (A similar analysis, but evaluating the derivatives at x = 1/2, 
yields p = 0.42.) 

The above analysis is complemented with a study of Pa,niCc), as suggested by Eq. (^. 
The graph of Pa,n{Cc) versus A„ shows (on log scales, see Fig. 5), a fair amount of curvature, 
making determination of f3 more difficult in this case. Linear fits to the data for n = 7-9, 8- 
10, and 9-11 yield, respectively, P = 0.471, 0.460, and 0.448, consistent with an approach to 
the value of 0.41 for large n. Verification of convergence must naturally await the evaluation 
of approximations for larger clusters. 

A further point of interest is the validity of the scaling hypothesis, Eq. (|^). The data 
collapse shown in Fig. 6, a plot of p* = A~^pa,n versus x = (C — Cc,n)/A„ provides support 
for the hypothesis. (In Ref. a similar collapse is demonstrated for the Domany-Kinzel 
model.) 



B. CAM analysis for i/n 



As shown in Refs. [T^-|T9[, the CAM approach is readily extended to dynamics. Let 



7.(0 = (7) 
p dt 

be the relaxation rate in the n-site approximation, and let 7(C) be the true relaxation rate. 
Then we expect 7 ~ |C ~ CcT"? while 7n ~ |C ~ Cc.nT"'*^^ in the n-site approximation, where 
the mean-field exponent is i^\\,aif = 1 for models with an absorbing-state phase transition 
| 27| . A scaling hypothesis, analogous to Eq. (|^), for the relaxation rate, is 




7n(c) = A::"^ ^ , (8) 
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where the scahng function g vanishes when its argument is zero. Supposing that g{x) ~ 
Ixl'^ii'*^^, we see that 7n(C) ~ 7nlC~Cc,n|, where the amphtude follows 

In ~ . (9) 

I determine the relaxation rate numerically for C, ~ Cc.n, and from these data extract the 
amplitudes 7„. The results, shown in Fig. 7, display substantial curvature on a log- log plot, 
so that direct determination of the critical exponent v\\ is not feasible. Simulations |TB|J2^ 
yield estimates for v\\ in the range 2.3 - 2.6. The CAM results are consistent with values in 
this range, if we include a correction to scaling. The solid line in Fig. 7 is given by 

ln7„ = 1.51nA„ + AA„-5, (10) 

with fit parameters A = 7.435 and 5 = 2.125, consistent with a correction to scaling expres- 
sion 7„ oc An" (1+AA„), with z/|| = 2.5. While results for larger clusters will be needed to 
determine v\\ with precision, one can at least assert that the present results are consistent 
with the rather imprecise estimates from simulations. 



VI. DISCUSSION 

I have devised a computational algorithm for generating n-site cluster approximations 
for a one- dimensional stochastic sandpile model with a fixed particle density. To facilitate 
the analysis I impose the height restriction zi <2. Analyzing the results for n < 11, I obtain 



the estimates Cc = 0.930(1), (3 = 0.41(1) and z/|| ~ 2.5, all in agreement with simulation p6 
While the results are not very precise, they provide significant independent support for the 
simulational findings, showing that FES models with a strict activity threshold belong to a 
universality class distinct from that of directed percolation (DP) [l2ll-p3[|. Sandpile models 



in which sites with an above-threshold height can remain stable (so-called "sticky grains"), 
have recently been shown to belong to the DP class |^2[, but such is not the case for the 



model studied here. Application of the methods used in this work to other sandpile models 
should prove illuminating. 

Acknowledgments 

I thank Mario de Oliveira, Tania Tome, Attila Szolnoki, Gyorgy Szabo, Geza Odor, and 
Paulo Alfredo Gongalves Penido for helpful comments. This work was supported by CNPq, 
and CAPES, Brazil. 



^ e-mail: dickman@fisica.ufmg.br 



8 



REFERENCES 



[1] P. Bak, C. Tang and K. Wiesenfeld, Phys. Rev. Lett. 59, 381 (1987); Phys. Rev. A 38, 
364 (1988). 

D. Dhar, Physica A 263 (1999) 4, and references therein. 

R. Dickman, M. A. Munoz, A. Vespignani, and S. Zapperi, Braz. J. Phys. 30, 27 (2000). 
e-print: |cond-mat / 99 1 0454 . 



M. A. Munoz, R. Dickman, R. Pastor-Satorras, A. Vespignani, and S. Zapperi, in Mod- 
eling Complex Systems, Proceedings of the 6th Granada Seminar on Computational 
J. Marro and P. L. Garrido, eds., AIP Conference Proceedings v. 574 (2001); e-print: 
|cond-mat/0011442| . 

R. Dickman, Physica A, in press; e-print: |cond-mat / 1 lOOiSj 

G. Grinstein, in Scale Invariance, Interfaces and Nonequilihrium Dynamics, NATO 
Advanced Study Institute, Series B: Physics, vol. 344, A. McKane et al., Eds. (Plenum, 
New York, 1995). 

C. Tang and P. Bak, Phys. Rev. Lett. 60, 2347 (1988). 

A. Vespignani and S. Zapperi, Phys. Rev. Lett. 78, 4793 (1997); Phys. Rev. E 57, 6345 

(1998) . 

R. Dickman, A. Vespignani and S. Zapperi, Phys. Rev. E 57, 5095 (1998). 

A. Vespignani, R. Dickman, M. A. Munoz, and Stefano Zapperi, Phys. Rev. Lett. 81, 

5676 (1998). 

M. A. Munoz, R. Dickman, A. Vespignani, and Stefano Zapperi, Phys. Rev. E 59, 6175 

(1999) . 

A. Chessa, E. Marinari and A. Vespignani, Phys. Rev. Lett. 80, 4217 (1998). 
A. Montakhab and J. M. Carlson, Phys. Rev. E 58, 5608 (1998). 

V. B. Priezzhev, J. Stat. Phys. 74, 955 (1994); E. V. Ivashkevich, J. Phys. A 27, 3643 
(1994); E. V. Ivashkevich,D.V.Ktitarev and V. B. Priezzhev, Physica A 209, 347 (1994). 
A. Vespignani, R. Dickman, M. A. Munoz, and S. Zapperi, Phys Rev. E62, 4564 (2000). 
R. Dickman, M. Alava, M. A. Munoz, J. Peltola, A. Vespignani, and S. Zapperi, Phys 
Rev. E64, 056104 (2001). 

M. Suzuki and M. Katori, J. Phys. Soc. Jpn. 55, 1 (1986). 

M. Suzuki, M. Katori, and X. Hu, J. Phys. Soc. Jpn. 56, 3092 (1987). 

N. Konno and M. Katori, J. Phys. Soc. Jpn. 59, 1581 (1990). 

A. L. C. Ferreira and S. K. Mendiratta, J. Phys. A26, L145 (1993). 

M. Rossi, R. Pastor-Satorras, and A. Vespignani, Phys. Rev. Lett. 85, 1803 (2000). 

R. Pastor-Satorras and A. Vespignani, Phys. Rev. E62, 5875 (2000). 

S. Liibeck, Phys. Rev. E64, 016123 (2001); Eur. Phys. J. B26, 75 (2002). 

S. S. Manna, J. Phys. A 24, L363 (1991). 

S. S. Manna, J. Stat. Phys. 59, 509 (1990). 

R. Dickman, T. Tome, and M. J. de Oliveira, e-print: |cond-mat/ 0203565 . 
J. Marro and R. Dickman Nonequilihrium Phase Transitions in Lattice Models (Cam- 
bridge University Press, Cambridge, 1999). 

W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical Recipes 
(Cambridge University Press, Cambridge, 1996). 

M. E. Fisher, in Proceedings of the International Summer School 'Enrico Fermi', Course 



9 



LI, (Academic Press, New York, 1971); M. E. Fisher and M. N. Barber, Phys. Rev. Lett. 

28, 1516 (1972). 
[30] T. Tome and M. J. de Oliveira, Physica A260, 99 (1998). 
[31] G. Grinstein, Z. W. Lai, and D. A. Browne, Phys. Rev. A 40, 4820 (1989). 
[32] P. K. Mohanty and D. Dhar, preprint: |cond-mat/ 0202345 . 



10 



TABLES 



iransition 


Probabiiity 


UzU lUl 


1/2 


— ^ zUU 


1/4 


Izu — > 2L)1 


1/2 


^ 102 


1/4 


^ 210 


1/4 


121 202 


1/2 


^ 211 


1/4 


220 ^ 202 


1/4 


^ 211 


1/2 


122 212 


3/4 



Transition probabilities for the restricted-height sandpile. Probabilitites are symmetric un- 
der reflection. 



n 




1 


0.5 


2 


0.75 


3 


0.80854 


4 


0.83682 


5 


0.85305 


6 


0.86378 


7 


0.87148 


8 


0.87736 


9 


0.88207 


10 


0.88594 


11 


0.88918 



Table II. Critical densities in the n-site approximation. 
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FIGURE CAPTIONS 

FIG. 1. Active-site density p versus D = maxc \dP{C)/dt\, for n=10, C = 0.8860. 

FIG. 2. Stationary active-site density p versus particle density (. Solid curves: n-site 
approximations for n = 3 - 11; points: simulation results for a system of 5000 sites. 

FIG. 3. Gritical point shift A^-^ versus 1/n. The solid line is a fit including a correction term 
as described in the text. 

FIG. 4. p'^ = dpa,n/dC, versus A„. The derivative is evaluated at x — l/A. The slope of the 
solid hue is -0.592, corresponding to (3 — 0.408. 

FIG. 5. Active-site density Pa,n(Cc) versus 5„. The slope of the solid line is 0.408. 

FIG. 6. Scaled active-site density p* — /S.~^pa,n versus scaled particle density x — {C,— 
Cc,n)/An, for n = 6 - 11. 

FIG. 7. Relaxation rate amplitude 7„ versus A„. The solid curve is a fit including a 
correction to scaling term as described in the text. 
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